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Abstract 

The incompressible Navier-Stokes equations currently represent the primary model 
for describing stratified turbulent fluid flows at low Mach number. The validity of 
the incompressible assumption, however, has so far only been rigorously established 
for adiabatic motions. Here, we show from first principl es that the use of available 
energetics and thermodynamics considerations applied to a turbulent mixing event 
associated with stratified shear flow instability r efutes the widespread idea that the 
incompressible assumption is also valid when diabatic irreversible effects are important. 
The main consequence is that dynamics and thermodynamics are strongly coupled in 
stratified turbulence. This departs strongly from the currently accepted wisdom, and 
calls for a complete revisiting of the physical processes governing stratified turbulence 
at low Mach numbers. 

The classical incompressible Navier-Stokes equations (INSE), often in conjunction with 
the Boussinesq approximation |Boussinesq( 1 903j| , currently represent the primary model for 
the description and understanding of fluid flows at low Mach numbers M a = U /c s , where U 
is the characteristic fluid velocity, and c s the isentropic speed of sound. However, despite its 
widespread use, the incompressibility assumption has been so far justified rigorously only in 
the context of adiabatic motions [Batchelor(1967) , Landau and Lifshitz(1987) Majda(1984) 



Lions and Masmoudi(199"8)] , the generalization to flows with diabatic and irreversible effects 



remaining an outstanding challenge Ansumali et al.(2005)Ansumali, Karlin, and Ottinger 



Physically, the primary motivation for the incompressibility assumption is arguably the sim- 
plification it brings about by decoupling the dynamics of the fluid from its thermodynamics. 
While this appears of interest for the study of adiabatic motions, for which thermodynamic 
effects are indeed of little importance, the advantage of doing so otherwise is less clear, since 
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it severely complicates the understanding of two important fundamental issues, namely: 1) 
whether the relaxation toward thermodynamic equilibrium plays any role in stratified tur- 
bulent flows? 2) whether the work of expansion/contraction B (see [TU] for a definition) may 
sometimes contribute significantly to the production of mechanical energy? Until now, these 
two issues have never been considered to be relevant to understand the behavior of strati- 
fied turbulence, and are therefore usually not discussed nor even mentioned in most current 
reviews of stratified turbulence, e.g., Gregg(1987), Fernando(1991)| . Arguably, this is so 



because until now fluid dynamicists have lacked a strong physical reason to question the va- 
lidity of the INSE to describe the behavior of turbulent fluid flows at low Mach numbers. In 
this letter, we present what we believe is the first strong case ever made against the validity 
of the INSE when irreversible diabatic effects are present, by showing that the assumptions 
underlying the INSE are inconsistent with exact energetics considerations based on the full 
compressible Navier-Stokes equations (CNSE). 

We start by pointing out the fact that, even though this is not universally accepted, the 
INSE fail to describe the laminar evolution of a thermally stratified fluid evolving under the 
action of molecular diffusion alone, both qualitatively and quantitatively. The main reason 
that the incompressible assumption does not make sense in that case is because the only kind 
of motion is due to the diabatic contraction/expansion of the fluid, which makes its velocity 
entirely divergent, with no solenoidal component. Moreover, it is also clear that the leading 
order energetics is inconsistent with that of the INSE, as internal energy (IE) is in that case 
the primary source of mechanical energy, with IE being first converted into kinetic energy 
(KE), and then ultimately into gravitational potential energy (GPE). Physically, the ultimate 
cause for the motion here is the departure of the fluid from thermodynamic equilibrium (TE 
thereafter), which molecular diffusion strives to restore, in accord with the prediction of the 
second law of thermodynamics stating that the entropy of an isolated system must increase 
until reaching its state of maximum entropy. 

If one agrees with this view, and hence with the primary role played by the work of 
expansion/contraction B in the laminar evolution of a stratified fluid, then it is natural to 
ask what could be so special about turbulence that could relegate the role of B to second 
order? Indeed, it is well accepted that entropy production (EP) considerably increases in a 
turbulent fluid, and it is even conjectured that EP could possibly reach its maximum possi- 



ble value permitted by circumstances Ozawa et al.(2001)Ozawa, Shimokawa, and Sakuma 



which is the well known principle of maximum entropy production (MEP). An increase in 
entropy production, however, suggests that a turbulent fluid evolves more rapidly toward 
its state of maximum entropy, i.e., toward thermodynamic equilibrium. Such a conclusion 
seems unavoidable, unless one can come up with some physical reason for believing that tur- 
bulence is able to generate some kind of counter-processes able to halt entropy production 
at some point before a state of maximum entropy is reached. We are skeptical that such 
a reason can be physically derived, given that EP can only stops when TE is reached, by 
which time the fluid would have become obviously nonturbulent. This idea departs from 
the currently accepted wisdom, which tends to consider that turbulent mixing homogenizes 
adiabatiatically conserved quantities. In the oceanic case, for instance, existing turbulent 
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mixing parameterization strive to homogenize potential temperature and salinity, whereas 
TE considerations would suggest homogenizing in-situ temperature and chemical potential 
|Fofonoff(1962)] instead. In any case, since the relaxation toward TE in the laminar case 
occurs in conjunction with the ultimate conversion of IE into GPE via the work of B, it 
seems legitimate to wonder whether in the turbulent case a higher EP rate could possibly 
mean a faster IE/GPE conversion via a correspondingly increased Bl Physically, this is not 
inconsistent with the idea that a turbulent nonlinear cascade toward smaller and smaller 
scales should increase the curvature of the in-situ temperature field T, which in turn should 
increase the magnitude of the local values of diabatic heating, with a corresponding increas- 
ing effect on the local velocity divergence and ultimately on B itself. In the classical view 
based on the INSE, this scenario is eliminated by assumption, but as said above it has never 
been discussed nor disproven. In fact, the main suprising result of this letter is precisely to 
suggest that this scenario is precisely what characterizes the behaviour of stratified turbulent 
flows. In order to establish our results, we take as our starting point the fully compressible 
NSE written under the form: 

p^- + VP = —pgz + z/V 2 u (1) 



Dt 



^ + P V-u = (2) 
DT, Q 



Dt T 

dl dl 



(3) 



r = iV,v), T = — , P = -- (4) 

where u = (u, v,w) is the three-dimensional velocity field, P is the pressure, p is the fluid 
density, v = 1/p is the specific volume, g is the acceleration of gravity, £ is the specific 
entropy, v is the molecular viscosity, I is the specific internal energy, T is the absolute 
temperature, while the local diabatic heating is given by 

Q = --V-F 9 + £ (5) 

where F q = —pC p nVT is the molecular diffusive flux of heat, with C p being the specific 
heat capacity at constant pressure, and k the molecular heat diffusivity, whereas e is the 
local dissipation rate of kinetic energy. Note that in writing ([1]), the Stokes' hypothesis of 
vanishing bulk viscosity is made R.E.Graves and Argrow(1999)| for simplicity, but this could 



be relaxed by modifying appropriately the definition of the dissipation of kinetic energy in 
the following. 

For simplicity, the fluid is assumed to be thermally and mechanically isolated (except 
for the work of the surface atmospheric pressure P a against changes in fluid volume). In 
that case, the evolution equations for the kinetic energy (KE) , gravitational potential energy 
(GPE), and internal energy (IE) reduce to: 

«m = - w+B - D - P «£ ( 6 ) 

dt dt 
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= w (7) 

dt v 1 

d(IE) 

where V Q \ is the fluid volume, while KE, GPE, and IE refer to the volume-integrated kinetic, 
gravitational potential, and internal energy respectively, 

u 2 



KE= / p—dV, GPE= / pgzdV, 
Jv 2 Jv 



IE= / pldV (9) 
Jv 

whereas the energy conversion terms W, B, and D are given by: 

Dv 



W 



J v P 9 wdV } B = J v P^pdV J 



D= / pedV. (10) 



v 



The focus here is on the energetics of an unstable parallel stratified shear flow of the kind 
that has been extensively studied by means of laboratory and numerical experiments, e.g., 
K.G.Winters et al.(1995)K.G.Winters, P.N.Lombard, J.J.Riley, and d'Asaroj|Peltier and Caulfield(2003)| 



and references therein. A typical evolution involves three stages. First, a laminar evolution 
where the shear flow and stratification are eroded through molecular viscosity and diffu- 
sion. Second, an unstable evolution associated with the development of turbulent dissipative 
structures and intense mixing. Last, a return to laminar evolution once the conditions for 
instability have been removed. Insight into the energetics of such a mixing event can be 
obtained by integrating Eqs ([ME]) over the relevant time interval, which yields the following 
budget equations for KE, GPE, and IE respectively: 

ARE = -W +B -D - P a AV ol (11) 

AGPE = W (12) 

AIE = D-B (13) 

where the overbar and A(.) represent respectively a time integral and difference over and 
between the endpoints of the time interval considered. Thus, W, B, and D are the to- 
tal density flux, work of expansion/contraction, and KE dissipation rate. It is custom- 
ary to measure the ratio ^mixing = AGPE/D, known as the "mixing efficiency". From 
TT2l) . one has W = ^mixingD, which is a popular way to parameterize W Gregg(1987)| , for 



^mixing appears to be reproducible observationally, with a value of 0.2 often cited and used 
Peltier and Caulfield(20d3)j . In the classical incompressible view of stratified turbulence, 
it is generally accepted to regard the works of expansion/contraction due to B and P a AV a i 
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as small or even negligible compared to the density flux W or dissipation D. If so, the 
predictions of the incompressible model for the remaining terms of the energy budget are: 

AKE , . AIE , N 

___ ~ _(i + lmixing ), « 1, (14) 

AKE 



' (X ^ 1 mixing)- (15) 

In the following, we show that the fully compressible NSE yield distinct predictions for ( fT4l) . 
but similar for (fl5l) . It follows that accurate measurements of AKE, AIE, and D should in 
principle allow Eq. (fT4l) to provide the basis for an observational test of the incompressibility 
assumption, but this remains to be carried out. 

In order to make progress, we need an independent way to estimate the magnitude of the 
work of expansion/contraction B\ as regards to the work of atmospheric pressure against vol- 
ume changes, we don't question the validity of assuming it to be negligible. In this letter, ad- 
ditional information about B is derived by computing separate budgets for the available and 
un-available components of GPE and IE, as defined by Lorenz [Lorenz(1955)]. Such an idea 



was originally used in the present context by K.G. Winters et al.(1995)K.G. Winters, P.N.Lombard, J. J. Rile 



in the framework of the Boussinesq approximation. Since we question the validity of the lat- 



ter, the work of K.G. Winters et al.(1995)K.G. Winters, P.N.Lombard, J.J.Riley, and d'Asaro 



is generalized to the fully compressible NSE. To that end, let us recall that the concept of APE 
is defined as the difference between the PE of the actual state minus that of a reference state 
which is defined as the state minimizing the total potential energy PE=GPE+IE in an adia- 
batic re-organization of the fluid parcels [Lorenz(1955) , K.GWinters et al.(1995)K.G. Winters, P.N.Lombar 



In mathematical terms, the reference state is defined by an invertible mapping taking a par- 
cel at the position x to its position x r = x r (x) in the reference state. By definition, such a 
mapping conserves the mass and entropy of the fluid parcels, which implies: 

S(x r .,t) = S(x, t) = T, r (z r ,t), (16) 

p(x r , t)J r = p(x, t) = p r (z r , t), (17) 

where J r = <9(x r )/<9(x) is the Jacobian of the transformtion. The reference state has several 
important properties which are straightforward to establish: 1) it depends upon z T only, 
which accounts for the last equality in (|T6|) and (TTT|) ; 2) it is in hydrostatic equilibrium at 
all times, i.e., dP r /dz r = —p r g\ 3) the velocity u r = (dx r /dt, dy r /dt, dz r /dt) of the parcels in 
the reference state satisfies the mass conservation equation Dv r /Dt = v r V r ■ u r . The main 
objective of available energetics is to write GPE and IE as 

GPE = AGPE + GPE r , IE = AIE + IE r 

i.e., as the sum of their available components (AGPE and AIE) and un-available (or reference) 
components (GPE r and IE r ), and to derive individual evolution equations for each, as shown 
below. With regard to GPE r , its definition is 



GPE r = / p r gz r dV r = / pgz r (x.,t) dV, 

JV r JV 
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so that its evolution equation is given by: 

d{GPE r ) 
Jt " 



J^pgw r dV = W r (19) 



by using the definition Dz r /Dt = w r and the assumption of mass conservation of a fluid 
parcel D(pdV) / Dt = 0. With regard to IE r , its definition is 

IE r = f p r I(Z r ,v r )dV r = f pI{H,J r v)dV (20) 

In order to derive an expression for the latter, we use the differential expression dIE r = 
T r dYj r — P r dv r = T r dT, — P r dv r . As a result, it can be shown that: 

d -^>=D-H-B r , (21) 

where 

H = j v p(l-^ QdV, B r = jf Pr^pdV, (22) 

by using the result that by assumption, DT, r /Dt = DT,/Dt = Q/T. Physically, we interpret 
B r and H as the fraction of the internal energy that can be converted into un-available and 
available GPE respectively. In the case of H, this is motivated by this term resembling the 
classical Carnot formula where the coefficient (1 — T r /T) plays the role of the classical Carnot 
efficiency factor, and pQ the role of the heating source. In the case of B r , an important result 
is that it can be rewritten as 

B ' - Du 



= f Pr-^-prdVr = - f U r ■ V P r dV r = W r 

JVr Dt JVr 



by first rewriting B r as an integral in the reference space, and then using integration by 
parts and accounting for the abovementioned three properties of the reference state. Once 
the evolution equations for IE r . and GPE r are known, it is straightforward to derive evolution 
equations for AGPE and AIE simply by subtracting those for GPE and IE minus those for 
GPE r and IE r , viz., 

fra.„„ «m =H _ B+Wrt (23) 

taking into account the above result B r = W r . We now return to the energetics of the mixing 
event considered above, by deriving separate budget equations for AGPE, AIE, GPE r , IE r 
as follows: 

AAGPE = = W - W r , (24) 
AAIE = = H-B + W r , (25) 
AGPE r = W r , (26) 
AIE r = D-H-W r , (27) 
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which in turn implies the following equalities: 



W = W r = B r = B - H (28) 

Eqs. ( 1251127]) and Eq. fl28|) are the central result of this paper, for they prove unambiguously 
that the "incompressible" view of stratified turbulence that the "compressible" work terms 
B aren H are small or even negligible compared to W is simply impossible. Indeed, the 
result that B — H = W in Eq. (|28[) imposes that either or both B and H be of comparable 
magnitude as W, which is sufficient to refute completely the validity of the incompressible 
assumption in presence of irreversible effects, regardless of the Mach number or of the gas 
or liquid nature of the fluid considered. On the other hand, we see no reason to question 
the validity and accuracy of the incompressible approximation applied to adiabatic motions. 
In fact, the latter can be applied to the estimation of H, since the difference between T r 
and T is entirely due to adiabatic compressibility effects, so that \T — T r \ = 0(T\P — P r \), 
where T = aT/(pC p ) is the adiabatic lapse rate (with a the thermal expansion, and C p 
the specific heat at constant pressure). In laboratory conditions at atmospheric pressure, 
T — T r is not expected to exceed a few mK for a fluid such as water or seawater for which 
r = O(10- 7 K/Pa) [Feistel(2003)1 . Although by no means a definitive proof, this suggests 
that the above relations can be simplified by neglecting H compared to B, in which case we 
expect B »s W. Moreover, the new predictions for Eqs. (|T4|) and f|T5|) are given by: 

AKE _ AIE _ 

~ 1) =0 ~ (1 ^mixing) j (29) 

AKE 1 ^ ( 

A Tip ~ — i _ ~ — (1 + Imixing + • • •) (30) 

J- 1 J- jmixing 

As stated previously, these new predictions for the ratios AKE/D and AIE/D (Eqs. 1291) 
differ substantially from those resulting from the Boussinesq approximation (Eq. UM . Indeed, 
such ratios would differ by about 20% for a mixing efficiency ^mixing = 0.2, the discrepancy 
between the two models increasing with the measured mixing efficiency On the other hand, 
both models are found to yield similar predictions for the ratio AKE/ AIE (Eq. [TBI versus 
Eq. (30D, so that only the ratios AKE/D and AIE/D would be useful to discriminate 
between the classical incompressible theory and the new compressible one presented here. 
Whether this can be done with present measurement capabilities is an issue that is beyond 
the scope of this paper, but that is would certainly be of interest to pursue. If not, an 
alternative would consist in comparing direct numerical simulations of a turbulent mixing 
event in the context of the compressible and incompressible NSE respectively. Such a project 
is currently underway. 

In summary, this letter refutes the validity of the incompressible assumption, which is oth- 
erwise well established for adiabatic motions, to describe fluid flows affected by irreversible 
diabatic effects. There does not appear to be a physical basis, therefore, for regarding dy- 
namics and thermodynamics as being decoupled in stratified turbulent fluid flows at low 
Mach numbers, in contrast with the currently accepted wisdom. The present results are 
important, because they call for a complete revisiting of the accepted ideas regarding how 
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stratified turbulence operates. Some of the several consequences implied by our results are 
the basis for |Tailleux(2007)| , and concern the so-called ocean heat engine controversy, and 
whether turbulence should be regarded as speeding up the convergence toward thermody- 
namic equilibrium, with important implications for turbulent mixing parameterizations in 
numerical ocean models. 
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